ومع‌صمتنو آهمتامصمط) ۱۲2 ولسبهه10 


۱ 1-9 0۴۰ ,2022 تمصحصی مه موی با بط ,2 ,۷۵ 


90 ز9 / 120]:10.22067 


به سوی علوم ریاضی ملمنن ومد 


۵ ک ۵ صصنا[خمها //۵9۰) 9( 


۳۸۳۸0 )501,۷۲۲ ۳۵۳ (۱۷۲۳۲/۱۳۱۵1: ۲۱۲۳۲۳۲۵۲۳۱۸۷۲ ۸۱ 
۵ ۳۱۳۰۵۳۳۱۲۳۷۲ 0۷ ۷۷۲۲۳۲ ۳۳۵۳۲۲,۲/۱۱۷۲۵ظ باص۳ ۷ ۲۳۱۷ 


51۳/1۳۲ 5۳1۳1۸۵۲۳۴۸2۲ ۸۳۲۱ ۳۱088۳] ۴ 


صق؟ رمعطه؟ واتممنصتا تلحععطمظ 0نطعطه عمصمنمگ لمعنمصمطنها که وتمو۲ ری‌تامصمطما که تصمصهممط 1 
0۵90۱۵ مهم جر تم امه .و 
اه دا وله هط 


جمص لصا مطه طمت صماطامتص مقعمصا متامطاه‌تهص مط رتممرهص فنطه ص. ماع جاوما 
۲ ۲36 مصتفنه مومع مط مامع ۲۵ م4صج 0معتن0 0ص مد تمامصصهتهمر معتنا۵و 
0۵ عمط مصتصصعهام مه وتحععم‌مظ فد ذ رصصممامعصر فتط معامع ۰10 ,0مط0فحط مصن)ناوره 
صملدامهر مط مصتمام‌هومی ها مه له مرصمتتصم مهن صه مصتفنه مامصصنهم 221 
مصتمد 5۵17۵0 ود صمامام‌جم‌حوانه هم رقصصماطاممطایاه مامرصته ومص از مصته هم 2۳0 
فص همحمت‌طه ۶ه ممارصعی آهه‌تمصصه مصتل‌زمم مانطاه رمعلظ .قحصمکعصه متتنوو۲ 
و حصماوامور مقعمص »تامجههم فطع متام عم 0مطتامص 0عفممرمهر مط رعع‌صهم۵ه۲ت هه 
910 


03 ,65]32 ,65712 .صمتای‌صتععهان هعمزطادو متامصمط) ]۷ 2020 
صتمتعصها متم۲ بل‌مطمظ مصت‌تاووه تماهووم مومع مععمص متامجاهتوظ .قل ۵ مک 
13-7-02 مصنلوی ماطماتوهظ 28-6-2022 ۸6۵۵۵۲6۵0 25-41-2022 ۲۵۷71۵۵0 4-12-2021 ۲۸۵۲۷۵0 ماع[ 
یرومم 7 

۹۱ 


۱ دوره ۰۲ شماره ۰۱ بهار و تابستان ۰۱۴۰۱ صص. ٩۱-۱۰٩‏ 
022.69. :1022067/00 .170۵1 


مقاله علمی-ترویجی 


به سوی علوم ریاضی 


۵ ۴ص -زقمتا //05۰ 


یک روش کارا برای حل مسائل معکوس سهموی با شرط اضافی 


سیمین شکرپز! و حسین آذری"" 


! دانشکده علوم ریاضی, دانشگاه شهید بهشتی» تهران» ایران 
۵90۵ سم جر تم امه .و 


که حا وله هط 


چکیده. در اين مقال» مسأله معکوس سهموی با پارامتر منبع مجهول معرفی شده و با استفاده از روش 
جداسازی عملگر به حل مسأله پرداخته‌ايم. برای حل اين مسأله. در ابتدا با استفاده از شرط اضافی به تعیین 
پارامتر منبع پرداخته و سپس با جداسازی مساله و تبدیل ان به زیرمسائل ساده تر. هر زیر مساله با استفاده 
از تبدیلات فوریه حل می‌شود. هم چنین ضمن ارائه مثال‌های عددی کارایی و همگرایی روش پیشنهادی برای 
حل مسأله معکوس سهموی نشان داده می‌شود. 


۱ معرفی 


مسائل معکوس از سال‌های دور مورد توجه بوده است. اولین مطالعات درباره‌ی مسائل معکوس و بدوضع 
به نیمه اول قرن بیستم برمی‌گردد که موضوعات آن ها به مسائلی در فیزیک (مسائل معکوس در نظریه پراکندگی 
ذرات)» ژئوفیزیک (مسائل معکوس در کاوش‌های الکتریکی» زلرله شناسی و نظریه پتانسیل) ستاره شناسی 
و دیگر علوم وابسته بود. در سال‌های اخیر حوزه کاربرد نظریه‌ی مسائل معکوس تقریبا به همه علومی‌که از 
فا کی تاقوا کب ره د اوق ایسی: 

امروزه تعداد زیادی از نشریات علمی به طور مستقیم و پا غیرمستقیم مسائل معکوس و بدوضع را بررسی 
می‌کنند. مسائل هدایت گرما؛ در دستة معادلات دیفرانسیل با مشتقات جزئی از نوع سهموی جای می‌گیرند. 


۰ ,651۷۲32 ,65712 .صمتای‌صتععمان #مزطانو متامصمط ]۳ 2020 
کلید واژگان. مسأله معکوس سهموی» روش جداسازی عملگر تبدیلات فوریه. 
تاریخ: دریافت ۱۴۰۰/۹/۱۳ بازنگری ۱۴۰۱/۲/۵ پذیرش ۱۴۰۱/۴/۷ انتشار بر خط ۱۴۰۱/۴/۲۲ 
"نویسنده مسئول 
نحوه ارجاع به اين مقاله: س. شکرپز» ح. آذری» یک روش کارا برای حل مسائل معکوس سهموی با شرط اضافی» به سوی 
علوم ریاضیء ۲ (۰)۱۴۰۱ شماره ۰۱ ۰٩۹۱-۱۰۹‏ 
۹ 


حل مسائل معکوس سهموی ۹۳ 
در مسائل مستقیم هدایت گرما. هندسه. شرایط مرزی, شرایط اولیه و خواص ترموفیزیکی معلوم می‌باشند و 
هدف محاسبه توزیع دما در داخل دامنه جواب می‌باشد. در مسائل معکوس هدایت گرما. یکی یا تعدادی از 
این اطلاعات نامعلوم بوده و هدف تخمین آنها با استفاده از اطلاعات اضافی دیگر مانند دماهای اندازه گیری 
شده در داخل ناحیه جواب میباشد. 
به دلیل اهمیت و کاربردهای فراوان مسائل معکوس, در این نوشته به بررسی و حل عددی مسأله معکوس 
سهموی مربوط به تعیین پارامتر منبع با شرط اضافی می‌پردازيم که به صورت زیر می‌باشد: 


(۱۰۱) و۱2٩‏ > (۶ره). و(اره) + رش (0)0 + مولا وا 
با شرط اولیه 

(۲۰۱) تک 0 ,(0)2 < (0 ,)0 

و شرایط مرزی 

(۳۰۱) 0 ,()وو > (۶ ,0/1 ,(9)6 > (0,۶) 


و قوظ افاقی به یکی از مرها 
(۴.۱) زا را 0 انا 


یا 


(۵.۱) 0 ون مار / 


مفروض است که در آن (7 > > 1,0 > * > 0: (2,۷)) < ٩2‏ و 0 < 7 توابع ولا رو جو 
علاوه برای برخی اعداد ثابت 0 < م۰ هسته (11)2 در رابطه زیر صدق می‌کند» 


(۱.ع) م > هه(ه) | / 


این نوع معادلات در مطالعات مربوط به فرآیند هدایت گرمایی» قضیه کنترل و نفوذ شیمیایی کا ربردهای زیادی 
دارند (مراجع [۷۰۶] را ببینید). 

معادله‌ی (۱۰۱) برای شرح یک فرایند انتقال گرما با یک پارامتر منبع استفاده می‌شود و معادله‌ی (۴۰۱) درجه‌ی 
حرارت در نقطه *2 در زمان ۶ را نمایش می‌دهد. هم چنین با استفاده از شرط اضافی (۵.۱) می‌توان انرژی 


مطلوب (7/)1 را روی دامنه تعیین کرد. 


1 س. شکرپز و ح. آذری 

مثال‌های متفاوت دیگری می‌توان یافت که شامل شرط انتگرالی اضافی می‌شوند. مثلا انتقال گرما در یک 
لوله نازک در انتشار ذره در پلاسما [۱۹۰۱۰]. برای مشاهده کا ربردهای بیشتر از این مسأله و مسائل مشابه 
مربوط به تعیین پارامتر منبع در فضاهای یک بعدی و ۱-بعدی به [۲۸-۲۴۰۱۵۰۸] مراجعه کنید. 

وجود. یکتایی و وابستگی جواب‌های این مساله در مراجع [۳۳۰۷۰۵] به اثبات رسیده است. 

روش‌های زیادی برای حل عددی اين مسأله وجود دارد» مانند تفاضل متناهی عنصر متناهی» عنصر مرزی» 
طیفی و حجم تتاهیه باء عتران فال فن [۵ | زوش خاض مشتاهی دراععحل انم مساله به کار فته اس انا 
در مورد تحلیل همگرایی و پایداری صحبتی به میان نیامده است. در مرجع 1 ٩‏ روش اویلر پسرو برای حل 
مسأله با شرایط اضافی (۴.۱) به کار رفته است. هم چنین روش هم‌مکانی سینک برای حل مسأله معکوس با 
شرایط مرزی غیرموضعی در مرج [۳۰] به کار رفته است. با استفاده از این روش» مسأله به یک دستگاه از 
معادلات اشگرالی رلعا دی ده و هگرانی روش اقرات م قرو 

در مراجع [۲۱۰۱۸۰۴] شرایط وجود و یکتایی برای جواب‌های مسأله معکوس با دو پارامتر مجهول مورد 
مطالعه قرار گرفته است. هم چنین مسائل مشابه در مراجع [۰۱۶,۳ ۰۲۰۰۱۷ 0۳۵۰۲۹ ۳۶] مورد مطالعه قرار 
گرفته است اما درمورد جواب‌های عددی با شرایط اضافی انتگرالی صحبتی به میان نیامده است. 

در این مقاله. از روش‌های جداسازی عملگر برای حل مسأله معکوس سهموی مربوط به تعیین پارامتر منبع 
با شرط اضافی استفاده می‌شود. در این روش. با جداسازی معادلات پیچیده و تبدیل آن به زیر مسائل ساده‌تر 
هر زیرمسأله با روش‌های انتگرال‌گیری مناسب حل می‌شود که با شرایط اولیه به هم مربوط هستند. سپس با 
ترکیب جواب‌های حاصل از زیر مسائل جواب مساله اصلی به دست می‌اید. 

اولین نتایج حاصل در مورد جداسازی عملگر» مربوط به روش جداسازی لی- تروترا (مرجع [۳۲] را ببینید) 
در سال ۱۹۵۰ می‌باشد. در سال‌های ۱۹۵۰-۱۹۶۰ پیشرفت‌های زیادی در زمینه این روش‌ها در روسیه به 
وود امد که می‌توان از نتایج فش امه تس ما رک نانک پدسامازشگی اشار تیوه اسان ترتع کار 
پژوهشی در این حوزه کا ربرد روش‌های لی* و استرانگ" در حل معادلات دیفرانسیل معمولی و جزئی می‌باشد 
(مراجع [۳۴۰۳۱۰۲۸۰۲] را ببینید). در این روش‌هاء با به‌کارگیری روش‌های مرتبه بالاتر (مراجع [۲۲۰۱۴] را 
ببینید) خطای جداسازی کاهش می‌یابد. هم چنین با انتخاب مناسب تر تابع تکرار اولیه. روش تکرار بهبودیافته 
حاصل می‌شود که می‌توان مرتبه دقت را افزايش داد (مرجع [۱۳] را ببینید). 

این مقاله شامل بخش های زیر می‌باشد: در بخش دوم پیش نیازهایی را که در طول این مقاله به ان ها نیاز 
داریم بیان می‌کنیم. در بخش سوم روش جداسازی عملگر برای حل مسأله معکوس ارائه می‌شود. در بخش 
چهارم نیز با ارائه مثال‌های عددی کارایی الگوریتم پیشنهادی را برای حل مسأله فوق نشان می دهیم. 


ام مزر[ 


۱/۹ 


آمتصمص 


"موه 
هموزر[ 


مهو 


حل مسائل معکوس سهموی ۹۵ 
۲ پیش نیازها 


به منظور تحلیل خطای روش جداسازی عملگر برای حل مسأله (۱۰۱)-(۴۰۱) یا (۱۰۱)-(۳۰۱) با شرط 
اضافی (۰)۵۰۱ نظریه نیم گروه های 0/0 مورد استفاده قرار می‌گیرد. با به‌کارگیری تبدیلات 


,(6)دو9 ۴ ,2 5 91 > (ر)2 ,6 و2 0 (1)2,۷ > (2,۷) 
۳ 


(۱.۰۲ ۱7 ع (2) ,(9)2,1 + رو(1) + ولا << ولا 
۳ ح 1 ِ 0 ,)10 5 (2,0/ 


با شرط اضافی به صورت زیر 


۳,۸0 - ۷, > (۳۰۲ 

یا 

/ 10((0-( (۳ 

که ور ان 
(0 ,)2 - (ه) من > (ه) ور (#,ع)ع(0)م + (بعیه - )۶ - 9۵ 
۰( ,)2 () > / - ()۳ <- (۲) ر(2)2*,1 - (۶) <- (0۶) 


سپس نظریه نیم‌گروه برای مسأله (۱.۲) با شرایط مرزی همگن به کار برده می‌شود و جواب‌ها با شرایط مربوطه 


با فرض این که (,۷)2 جوابی از مسأله (۱.۲)-(۲۰۲) باشد» می‌توان نوشت؛ 
9( و(ق) له (6 ,)و + (۶,ه)یا(گ)م + (۲ر )مولا 
همچنین با درنظر گرفتن شرط انتگرال 


(۳.۲( ۰ > > . ,() ۷( «9)2,۷((4 + (2۵ + ویین)(11)2 / 


۹۶ س. شکرپز و ح. آذری 
لذا داریم 


(.۵) 7 > ۳ 0 و 5 تست تک )۳ ۳ ()2 


و به‌طور مشاپه برای شرط اضافی به فرم (۳۰۲) می‌توان نوشت: 


۹ ۳ - ,ول( )1 ۳ 2۳ 


لا ۳ .« 77 < ()2 
فرض کنیم [(6 ,)1 جد ه] < (0) ۷ و [(2,0)و + «] < (0)0. مسأله مقدار اولیه زیر را در نظر 
می‌گيريم؛ 
(۷.۲) بو  )2(‏ (0,۷ + ۸۷ ع< م۷ 
با > 0 و < (۲)0 


که در آن ۸ یک عملگر خطی بی‌کران است (در این جا ۸ مشتق دوم در جهت 7 می‌باشد) و 
,۲ بت 2 ۷« (0,1] : (۲ ,)۳ 
,(6)0 + ۲()م < (۷ ,)۳ 


تابع غیرخطی از متغیرهای / و (1) ۷ می‌باشد و ()2 با رابطه (۵.۲) يا (۶۰۲) ب‌دست می‌آید. در این رابطه» 
فضای جواب مسأله (۱۰۲) می‌باشد. 


۲ روش جداسازی عملگر در این بخش مسأله مقدار اولیه (۷۰۲) را روی بازه [0,1] در نظر می‌گيريم. 
فرض کنیم (1۳) ۲ تج ,1 جواب عددی در زمان 77 < " و با طول گام 7 باشد. در این صورت با به‌کارگیری 
روش جداسازی لی» زیرمسائل زیر را روی هر زیربازه [1۳۲۲,"/] در نظر می‌گيريم, 


ن 00 

(۸.۲) ک زر ۶ ,40 ثٍِِ 
0 
تا ۳ 8 < سب 
( هت( ,ماهس 

۳ 

9 سل 

13 5  )ر و < ‌ ۳ و(‎ )٩.۲( 
۰ )نس > ()ن‎ ( 


حل مسائل معکوس سهموی ۹۷ 
زیرمسائل فوق می‌توانند به صورت کارا حل شوند که در آن ,۷ تقریب جداسازی در ۷۳ < ۸ می‌باشد. با 
به‌کارگیری الگوریتم های عددی برای حل زیرمسائل .)٩.۲(‏ جواب می‌تواند به صورت زیر به‌دست بیاید. 

ویر زر ۲ ,( )سمل > (۳۲) ره 


,( )نله < (۷)0۳ ,() سره < ( )0 


که در آن زیر مساتل جداسازی با شرایط اولیه به هم مربوط می‌شوند. بنابراین جواب تقریبی به صورت 
)+۵ > 3 می‌باشد. 


۳ _ روش جداسازی عملگر برای حل مسأله منبع معکوس سهموی 


در این بخش روش جداسازی عملگر ۳ و محاسیه حواب‌های عددی در هر گام زمانی به کار می‌رود. 
با فرض مقادیر تقریبی (1,)م»04 روش جداسازی عملگر در 7-امین گام مقادیر (6) با استفاده از شرایط 
اضافی به صورت زیر محاسبه می‌شوند: 

2 9 تک 2 و تایهزا 55 )۷ 


۱0 ت(0 ب 7 


)1( 


یا 


۵( ,)0 (2) ۲ ۳ - ,ول( )1 1 ۳( 
9 


که در آن (#,:2) 1 جواب مسأله متناظر با شرایط مرزی همگن می‌باشد. در نتیجه با استفاده از مقادیر محاسبه 


0 0. م < (2)0 
شده (0)1 در امین گام زمانی زیرمسائل )٩۰۲(‏ به صورت زیر تعریف می‌شوند 


(۳.۲( و و) مرول ۳ ر ۳ ,)نله كِِ ِ 1 ۰ یف ,(9):1 ‌" ۲( 5 ( بآ ۶ ول 


و 
(۳.۲( ره رز کر ی ,820 < ۸40 < پر 


که در آن 1 - ۷ 2,۰۰۰ ,1 ,0 < ۶ و 2۷ تعداد نقاط شبکه در بعد زمان است. زیرمسائل فوق برای محاسبه 


مقادیر تقریبی )1 ار در (1 ِ 0)- امین گام زمانی حل می‌شوند. 


۳ تبدیلات فوریه برای حل زیرمسائل جداسازی برای حل زیرمسائل (۲۰۳)-(۴۳) با شرایط مرزی 
متناوب. تبدیل فوریه سریع به کار برده می‌شود. هم چنین در موردی که شرایط مرزی دیریکله و نامتناوب 
باشند. روش تبدیل فوریه سینوسی استفاده می‌شود. تبدیل فوریه کسینوسی نیز می‌تواند برای شرایط مرزی 
نیومن و نامتناوب استفاده شود. 


۹۸ بزی: شکرپز و ح. آذری 
مسأله مقدار اولیه (۳۰۳)-(۴۰۳) روی دامنه [0,1] < [1 ,0] را در نظر بگیرید. با فرض این که شرایط 
مرزی متناوب باشند. تابع (0)2:,7 به صورت زیر نوشته می‌شود: 


0 


(۵.۳) در زان 


0-<و] 


که در آن ضرایب فوریه به صورت زیر هستند: 
(۶.۳ ۳ > (۵)1 
0 
فرض کنیم شبکه یکنواخت از 1 + 1 نقطه به صورت زیر را داشته باشیم 
<< ر حد < رل 


با شرایط مرزی تناوبی» داریم (۲ ,20 ‌- (۲ ,زرتد) نا بنابراین )1 35) ۸ در نقاط و می‌تواند به صورت 


زیر تقریب زده شود 


۷-1 
۳۳( (۵)+۵ ر << (8 روت)ن 
1-0 
که کاخ 
۷۲-1 1 
۸.۳) و( روت را جر - 6۵ 
(0<ز 


این رابطه تبدیل فوریه گسسته تابع (2,۲) نامیده می‌شود و می‌تواند برای محاسبه مقادیر تقریبی ( ,)0 
ه‌کار برده شود. 

مقدار ون را بدعتوان حواب اولیه معادلهی (۳:۳) در ومان قابت ۶۳ در نظر بگیرید که برایر مقدار اولیه ول 
می‌باشد. با به‌کارگیری تبدیل فوریه سریع؛ داریم 


۱ ۷-1 ۱ ۷-1 ۱ ۷۲-1 ۳ 
(۵.۳) 5 ت "1 *() برزه 3 (۷) ِِ ۶و ()یرزه 5 
0-<:/ 1-0 (اک 
و00 < (2)0 


و در نتیجه. 


۵ )( - )0(۵۸)۵ + 9)0( )۱۰۰۳( 


حل مسائل معکوس سهموی ۹۹ 
با حل معادله (۱۰۰۳) می‌توان نوشت.» 


#۷ 
(۱۱۳) ,1 ۸ , ...0,1 ول(و)ب(2(۹۵ولی / + پرونم ۲۵۹8 ملع - ()یر۵ه 
0 


که در آن ,420 ضرایب فوریه )1:0 هستند و با رابطه (۸۰۳) محاسبه می‌شوند و (8). ضرایب فوریه تابع 
(5 ,)9 می‌باشند. انتگرال‌ها در رابطه (۱۱۰۳) با استفاده از قاعده سیمپسون تقریب زده می‌شوند و جواب 
ای هراق که کشا دای مایق ماک بعش به نع 
(1 ,)نا 

به طور مشابه, برای محاسبه جواب معادله‌ی (۴.۳)) دا ریم: 


۷۲-1 
 2)۵( ۳, ۳۳۳‏ < ع (۶,ه)ن 
۲<0] 
و با حایگذاری در معادله‌ی (۰)۴۰۳ می‌توان نوشت. 
01 
۳۳ ,(72)1 ۳ - << 2 
00۰ 5 (0:)0 


با حل رابطه‌ی اخیر داريم 
(۱۴.۲ 1 0 ها رنه < (0)6 
بنابراین جواب (0)2,۶ محاسبه می‌شود که متناظر با مقدار (***/,) می‌باشد. با به‌کارگیری فرمول 
جداسازی عملگر لی؛ داریم: 
و مقادیر (1,:)مه" در (1 ۳ 1)-امین گام به صورت 
ی پر و 
محاسبه خواهد شد. 
برای حل مسأله منبع معکوس فرضی با شرایط مرزی دیریکله و نامتناوب 
وا وا 0 ها با که رما 
تبدیل فوریه سینوسی و برای شرایط مرزی نیومن 


تبدیل فوریه کسینوسی به کار می رود (مرجع [۲۳] را ببینید). 


۱۰۰ س- شکرپز و ح. آذری 
مسأله (۳۰۲) را روی 0,7 >«( < م۲) با شرایط اولیه و مرزی مفروض در نظر می‌گيريم. با استفاده از 


سری فوریه سینوسی. جواب می‌تواند به صورت زیر نوشته شود 


2) 0 51) ( )۱۵۰۳( 
۳۳ 

ی مروت یقرت 3 م وشوو 
1 

012 109 511) ( )۱۶۰۳( 
۳۳ 


و در آن ضرایب فوریه است و داریم ۷10 < (۰10)2:,0 در رابطه‌ی اخیر تبدیل فوریه سینوسی به صورت 


زیر خواهد بود» 
1 

۱۷۰۳ ۰( ۲۵۲ )تلو () ,)نب 2 (0)6 
0 


در نتیجه با جای‌گذاری (۱۶۰۳) در رابطه‌ی (۳۰۳) می‌توان نوشت. 


ِ 


۵ 510)72( - 0۲( ۲ ,(7۲2)ط51 2۵ + (172)و ( )رما‎  )۱۸۰۲( 
۳ ۲-1 
۵:)0( < 00۸, 
که به رابطه‌ی زیر منجر می‌شود.‎ 
رو < (۵:)0 ,(9)6 + () ,را( > (6) ماس‎ )۱٩۰۳( 
جواب معادله‌ی اخیر می‌تواند به صورت زیر نوشته شود‎ 


1 
(۲۰۰۳) و نت ,2 1 | ,05 (5) م2 2)0(۵ وم / ِ پر *2)5(9 وم تِ ()0۸ 
0 


که در آن 420۸ ضرایب فوریه تابع (100)2 هستند که با رابطه‌ی (۱۷۰۳) محاسبه می‌شوند. انتگرال‌ها در رابطه‌ی 
(۲۰۰۳) با استفاده از فرمول سیمپسون تقریب زده می‌شوند. لذا جواب معادله‌ی (۲۰۳) محاسبه می‌شود که یک 
مقدار اولیه برای (۴۰۳) می‌باشد یعنی (#,:) < (0 ,)۰۷ 

به طور مشابه. برای محاسبه‌ی جواب معادله‌ی (۴.۳)» داریم 


(۲۱۰۳) ب(167۲2) او ()بر0 ( 7 )0 


حل مسائل معکوس سهموی 1 
و ضریب فوریه ی (8:)6 با رابطه‌ی زیر به دست می آید» 
1 
(۲۲۰۳) ۰( )له (6 ر) 0 2 06 
0 


با جایگذاری (۲۱۰۳) در معادل‌ی (۴۰۳) داریی 


(۲۳۰۲) ,(720)1 ۶ - - ی 


که در آن. جواب به صورت زیر محاسبه می‌شود. 
(۲۴۰۳) و ,۳*۵ کی -- (0)۶ 


در نتیجه (6 ,0)2 محاسبه می‌شود که متناظر با ( 2,۳۲۴ می‌باشد. با استفاده از فرمول جداسازی لی» داریم 


و در نتیجه مقادیر (,2)مع در (1 + 0)-امین گام محاسبه می‌شوند. 


۴ نتایج عددی 


در این بخش برای ارزیابی صحت روش جداسازی عملگر به مرور چندین مثال عددی می پردازيم و هم چنین 
زیرمسائل با استفاده از تبدیل فوریه حل می‌شوند. روش تبدیل فوریه سینوسی برای بسط جواب معادلات 
(۳۰۳)-(۴.۳) به کار برده می‌شوند و جواب‌ها در فضای فوریه در (1 ۳ 11)-امین گام زمانی محاسبه می‌شوند. 
وان اثبات کارایی روش جداسازی فلگ لی: نتایج عددی را مقادیر مختلف 7 و محاسبه کرده آیم. 

فرض کنیم 26 فضای جواب مسأله (۱۰۱)-(۳۰۱) با شرایط اضافی مفروض به فرم (۴.۱) یا (۵.۱) باشد 
گنک هام باناع اک مریم |۱۷ زا سا رای شراب هی بت ترس فضای ردام عفر 
جواب‌ها به صورت زیر محاسبه می‌شود. 


(۱.۲) را( موت)موله - (1 ,)| < برااز6|| 
که در آن ره متناظر با بت ( > > 1) می‌باشد و (7 ,)مره جواب تقریبی محاسبه شده با استفاده از 


روش جداسازی در مسأله معکوس می‌باشد. 


در این قسمت نرم ماکزیمم خطا به صورت زیر تعریف می‌شود» 


(۲.۴( زا و یلاع( ۲ و ۳ < 0۳0 - گر[ 


۱۰۲ س. شکرپز و ح. آذری 


و مرتبه محاسباتی روش پیشنهادی عبارتست از 


هم چنین خطای پارامتر منبع در زمان انتهایی به صورت ((1) - ( )۱۸ < (6,)1 تعریف می‌شود که 
هو اراز )جوا کفرییی امشا: 
مثال ۱۲:۲: فسالة (۱۰۱)-(۴۰۱) روی دامنه‌ی [1 ,0 < [1 ,0 7( با داده های مفروض زير را در نظر 
می‌گیریم» 


)ده > (ه)ون 

2 

,۵ < (2)8و 6 2 ()9۱ 

1 2 ۶ 

9 ۳۱۸ .2 ۸ 
و جواب‌های دقیق عبارتند ان 
7 2 
,گم - 1 << () ,)005 ۳ 


در این مثال» خطای پارامتر منبع در زمان انتهایی به ازای مقادیر 0.01 27۳ (/ قتیاو فیشیرته اد 10 < 1۰0 
نرم ماکزیمم خطا و مرتبه محاسباتی روش جداسازی عملگر لی در جدول ۱ در زمان انتهایی 1 برای 0.01 < «/ 
و مقادیر مختلف 7 محاسبه شده‌اند. 


جدول ۰۱ نرم ماکزیمم خطا و مرتبه محاسباتی روش جداسازی لی در مثال ۱.۳ برای 
1 - «] و مقادیر مختلف 7 


7] روش تبدیل فوریه 
نرم ماکزیمم مرتبه محاسباتی 

0.3708 3 

10212 0.1765  |دج‎ 

12-133 0093 7 

27 0.0912 1-3409 
و 0.0098 10081 
مج 0.0025 1-334 
وج 0.0005 30 
| 0.0002 10092 


شکل ۱ جواب‌های تقریبی 2 )ون و تابع خطای مربوط به روش جداسازی عملگر را در دامنه‌ی : ) 1 


(1 > ۷ > 1,0 > 2 > 0 برای 0.01 < 7۱ و مد < 7 نشان می دهد. 


ناهن 6۲۱۵۲ 


(۳) 


شکل ۱ جواب‌های تقریبی ره اس )( و تابع خطای مربوطه (ب) در مثال ۱.۴ ترا 
1 2 ۸ و بط < 7 


6۲۲۵۲ ۲ 


104 103 10۶ 101 10 
0 


شکل ۲: نرم خطا در مثال ۱.۴ برای 0.01  <‏ و مقادیر مختلف 7 


نرم ماکزیمم خطا با استفاده از روش جداسازی عملگر لی در شکل ۲ و به ازای مقادیر مختلف طول گام 
زمانی نشان داده شده‌اند. در این شکل می‌توان دید که با کاهش طول گام زمانی» نرم خطا نیز کاهش می 
یابد که بیانگر همگرایی روش پیشنهادی می‌باشد. مقایسه نتایم حاصل با موارد مشابه در مراجع [۱۱] و [۱] 
بیانگر برتری روش فعلی می‌باشد. در این مراجع روش تحلیل هموتوپی و روش کرانک نیکلسون برای حل 
مسأله (۱۰۱)-(۴.۱) بکار برده شده است و کمینه مقدار نرم ماکسیمم خطا * 10 می‌باشد. هم چنین در مراجع 
ذکر شده نرخ همگرایی در مثال‌های عددی بررسی نشده است. هم چنین با توجه به بکارگیری روش جداسازی 
عملگر هر زیر مسأله با روش بهینه متناظر با عملگر موجود حل می‌شود و لذا زمان محاسباتی بکار رفته در 
روش پیشنهادی کمتر است. 


منال .۳ تیاه (۱۰۱)-(۴۰۱) را روی دامنه‌ی 1 00 ۴ [1 ,0 ت 0 با داده های زیر در نظر می‌گیریم» 


کم ((م )هه 2۵ - 0 + (ماصن (1 - 6)) > ۶2,0 


+ 511)2(6 ( 37 )وا 


۱۰۴ س- شکرپز و ح. آذری 


1 

,)هه <- < (0)و9 . ,20 ():9 
۳ 1 

ان بر < 9 


و جواب‌های دقیق عبارتند ان 
1 
.11 - 3- << (1) ع() دز و 03 
روش حداسازی لی برای حل این مثال استفاده می‌شود و با استفاده از روش تبدیل فوریه زیرمسائل جداسازی 


حل می‌شوند. در جدول ۲ ۰ مقادیر نرم ماکزیمم خطا برای جواب‌های (2:,۶),ء" در زمان انتهایی به ازای 
1 - ۲ و مقادیر مختلف 7 ارائه شده‌اند که پیان گر کارایی و توانایی روش جداسازی ی بط تولید 


حدول . نرم ماکزیمم خطا و مرتبه محاسباتی در مثال ۲.۴ برای 1 << (/ و مقادیر 


مختلف 7 
7 | روش تبدیل فوریه 
نرم ماکزیمم مرتبه محاسباتی 
3 0.2026 
ِ 10035 10893 
2 00495 1043 
7 0.0240 1.0430 
و 00113 13408 
0.0049 1-0387 
5چ| ‏ 0.0018 10294 
0.0005 1104 


می دهد. خطای پارامتر منبع در زمان انتهایی با مقادیر 0.01 < 2 و مش < 7 عبارتست از 4 10 ۶« 0.5. 

مقادیر نرم خطا برای جواب‌های (2,۶),,» در زمان انتهایی روی فضای باناخ 26 در شکل ۴ آورده شده 
است. همانطور که مشاهده می‌شود با افزایش تعداد نقاط شبکه مقادیر نرم خطا کاهش می یابد. مقایسه 
نتایم حاصل با موارد مشابه در مراجع [۱۱] و [۱] بیانگر برتری روش فعلی می‌باشد. در این مراجع روش 
تحلیل هموتوپی و روش کرانک نیکلسون برای حل مساله (۱۰۱)-(۳۱) بکار برده شده است و کمینه مقدار 
نرم ماکسیمم خطا 103 می‌باشد. همچنین در مراجع ذکر شده و۳ همگرایی در مثال‌های عددی بررسی نشده 


انشتاه 


منال ۳.۴ ی له (۱۰۱)-(۳۰۱) و (۵۰۱) را روی دامنه ,0 > [1 ,0 تک 7 و با داده‌های زیر در نظر 
بگیرید» 


و((2)صنه (1 - ه)( 2 +01 - (تاومه 2-)* 6 < (2,۶) 


)هه (1 - )> (ن)ونه 


سل بقل کوب بهیزی ۵ 


)1( (ب) 


شکل ۰.۳ جواب‌های تقریبی را )( و تابع خطای مربوطه (ب) در فا ۳ برای 
1 2 ۵ و بط < 7 


100 
10 
3 
ِ 
9 
4 
5 
0۳۹ 
نامز سهب 
: ۲ 104 
10 101 10 103 104 


0 


,0۰ < (۶)وو < ()9۱ 
«(1 - 1طلع)*7» < ()7] ۳ ۱۱۵ 
جواب‌های دقیق عبارتند ان 
,مب 1 < (۶) ر(27)صله (1 - 2)* 6 < ( ,2اه 

روش جداسازی لی برای 1 < 7 در اين مثال به کار برده می‌شود ونتایج محاسبه می‌شوند. شکل ۵ جواب 
تقریبی و تابع خطای مربوطه ر در زمان انتهایی به ازای 1 - ,| و سس << 7 تفای می دهد. این نتایج 
صحت بالای روش را براي حل این مسأله نشان می دهد. 

در این مثال» نرم خطا برای پارامتر منبع در زمان انتهایی 1 7 عبارتست ارگ 10 < 3.)09. به علاوه نرم 
های مختلف خطا برای جواب‌های (71 وا در زمان انتهایی روی فضای باناخ ۳۹ در شکل ۶ آورده کرهش : 


۱۶ مر : شکرپز و ح. آذری 


(ا,) 


(۳) 


شکل ۰۵ جواب‌های تقریبی (1,:)م. (۱) و تابع خطای مربوطه (ب) در مثال ۳۰۴ برای 
1 << و بط < 7 در زمان 1 < 7 


0۲۵۲۵۳0 
2 
0 


104 103 102 101 10 
‌ 


۳ نرم خطا در مثال ۲۰۳ برای 0.01 < / و مقادیر مختلف 7 


ین نتایج برای 0.01 -- :۸ و مقادیر مختلف ۲ محاسبه شد‌اند. در این شکل با افزایش تعداد گام های زمانی؛ 
نرم خطا کاهش می یابد که همگرایی روش را نشان می دهد. 
این مثال در مرجع شماره [۳۷] با استفاده از روش هم‌مکانی سینک حل شده است. در این مرجع با 
بکارگیری تبدیلات خاص معادله به معادله انتگرال ولترا تبدیل می‌شود و سپس با بکارگیری یک روش صریح 
و خصوصیات توابع سینک این معادله با یک سیستم از معادلات جبری خطی جایگزین می‌شود. مقدار نرم 
ماکسیمم خطا برای جواب تقریبی در مرجع [۳۷] مقدار * 10 ۲ 6 می‌باشد. در حالی که با استفاده از روش 
جداسازی عملگر مقدار متناظر عبارتست از *107 8 که بیانگر برتری روش پیشنهادی می‌باشد. 
۵ نتایج 
در اين مقاله روش جداسازی برای حل مسأله معکوس سهموی با شرط اضافی مورد مطالعه قرار گرفته 
است. برای حل این مسأله در ابتدا با استفاده از شرط اضافی به تعیین پارامتر منبع پرداخته و سپس با 
جداسازی مساله و تبدیل آن به زیرمسائل ساده تره هر زیر مساله با استفاده از تبدیلات فوریه حل می‌شود. 
درنتیجه با ترکیب جواب‌های عددی حاصل از زیرمسائل» جواب مسأله مفروض به دست می آید. نتایج حاصل 


حل مسائل معکوس سهموی ۳ 


از کاربرد روش جداسازی برای حل مسأّله منبع معکوس نشاندهنده کارایی و صحت روش پیشنهادی می‌باشد. 
هم چنین با کاهش طول گام زمانی خطا کاهش می یابد که همگرایی روش را نشان می دهد. 

در کارهای آینده می‌توان از دیگر روش‌های جداسازی مانند روش جداسازی تکراری برای حل مسائل معکوس 
استفاده کرد. هم چنین با در نظر گرفتن مسأله در ابعاد بالاتره می‌توان روش‌های جداسازی را برای حل آن 


بررسی کرد. 
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